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Abstract 

^■i Recent measurements of semileptonic i3-meson decays into r-leptons are somewhat higher than 

expected in the standard model. Although the deviations are less than 3a, they suggest the 



possibility of new physics affecting primarily the r-lepton. In this paper we examine these results 
within the context of non-universal left-right models. We find that strong constraints from b — > 57 
on W — W mixing lead to a prediction of approximately equal enhancements for the B — > Dtu and 
the B —7- D*tv modes. The model predicts approximately the same enhancement for the inclusive 
semileptonic rate B — > X c tv as well as for the leptonic decay B~ — > t~v. An enhancement in the 
leptonic decay B~ — > is also possible but is not uniquely correlated with the other modes. 
For this explanation to be viable, the mass of the non-universal W would be below 1 TeV, within 
LHC reach. 
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I. INTRODUCTION 



Recent measurements of B decay modes involving r-leptons have shown hints of devia- 
tions from the standard model (SM). In particular, for the semileptonic b — » crv modes, for 
which BaBar has reported recently that [l| 

, , B(B Dt-u t ) 
The corresponding numbers from Belle are ^ 



, N B(B Dt"v t ) 

R(D) = — V= — ^ = 0.35 ±0.11 

v ; B(B ^ D£-ue) 

, B(B D*t-v t ) 

Both the BaBar and Belle results are a bit high compared to the SM expectation (3), 0], 

R(D) = 0.297 ±0.017 

R(D*) = 0.252 ±0.003 (3) 

At the same time the discrepancy between theory and experiment for the leptonic B — )■ tv 
rate appears to be getting smaller, with the latest Belle result with hadronic tags being 
0.72ioi25±0.11) x 10~ 4 [sj]. Below we quote Belle's combination of this result with their result 
based on semileptonic tags [5j, as well as BaBar's number j^] and a recent SM prediction by 
the CKMfitter group 0- 

B(B -)• ru) BMe = (0.96 ± 0.26) x 10~ 4 
B{B Tu) BaBar = (1.79 ± 0.48) x 10" 4 

B{B t^sm = (0.719ig;#jj) x 10" 4 . (4) 

The apparent discrepancy in the semileptonic modes has inspired their examination 
within two Higgs doublet models. The data disfavors type II models 0] where it is not 
possible to simultaneously accommodate the enhancements in R(D) and R(D*) [l[. Type 
III 2HDM have been also studied and are consistent with the data [9J. 2HDM models are 



also discussed in this context in Ref. 10 . The results have also been discussed in the context 
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of R-parity violation [llj and more generally with the aim of making further predictions 
They have also inspired new calculations of the form factors (l3| . 

In this paper we consider another new physics possibility to address this enhancement 



In the context of non- universal left-right (LR) models [15J, these modes receive additional 
contributions from the right-handed charged currents that are suppressed or absent for 
modes involving the first two generation leptons. The model contains one light right-handed 
neutrino that would ultimately be responsible for the enhancements in the semileptonic 
5-meson decays by adding the channel ru R to the possible final states. 
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II. NON-UNIVERSAL LR MODEL 



The model we consider is a variation of LR models [l4j, and we have detailed it previously 



151 ] . Here we review the salient features relevant for our discussion. The gauge group is the 
usual SU(3)c x SU{2)l X SU(2)r x U{1)b-l- With this gauge group there are two charged 
gauge bosons, Wl and Wr for SU{2) L and SU(2) R , respectively. In general they mix with 
each other and the mass eigenstates W (light) and W (heavy) are a mixture of Wl and 
Wr. This mixing can be parameterized as 

W L = cos^wW - sin£ w W' , 

W R = sm( w W + coB( w W' . (5) 

The mass splitting between W and W' in this model is achieved by introducing the Higgs 
multiplets Hl '■ (1, 2, 1)(— 1) and Hr : (1, 1, 2)(— 1), and by letting their vacuum expectation 
values (vev) vl and Vr be different. In our notation for the quantum numbers of the 
multiplets, the three numbers in the first bracket are the representations under SU(3)c, 
SU(2)l and SU(2)r respectively. The number in the second bracket is for B — L. The 
Wl — Wr mixing can be induced by a bi-doublet scalar : (1,2, 2)(0) in which both of the 
diagonal fields in <fi have non-zero vev v\^- 

The non-universal version of the L-R model we consider here is the one proposed in 



Ref. (l5|. In this model, only the third generation interacts with Wr, and the first two gen- 
erations have SM like interactions. The first two generations (indicated by the superscript) 
thus transform under the gauge group as 

: (3, 2, l)(l/3) , U 1 / : (3, 1, l)(4/3) , D 1 / : (3, 1, l)(-2/3) , 
4> 2 :(1,2,1)(-1), 4' 2 :(l,l,l)(-2). (6) 

The third generation is different, as its right handed particles also transform under 
SU(2)r. The transformation of the third generation under the gauge group is given by 

Ql:(3,2,l)(l/3), Q|:(3,l,2)(l/3), 

L|:(l,2,l)(-1), L|:(l,l,2)(-1). (7) 

This model is anomaly free. It accommodates an enhanced interaction between third gener- 
ation fermions that becomes of electroweak strength for process that involve a fermion pair 
from the first two generations and one from the third generation. Its phenomenology has 
been studied in detail, with the most stringent bounds on the Z' arising from LEP2 [l5[ 
and most recently from the study of pp — > t + t~ at the LHC 1 (l6| . We are not aware of 
any LHC searches for non-universal W' bosons so far, but this paper would be an additional 
motivation to conduct them. 

In the mass eigenstate basis the quark-gauge-boson interactions are given by, 

C w = -^U L ^V KM D L (cosZ w W+ -sm£ w W'+) 

-^U R rV R DR(sint w W+ + cos£h/W/ + ) + h. c, (8) 



1 Analysis of this data within our model is not finished but we expect a constraint on the Z' mass around 
1 TeV, of the same order as those quoted by CMS for other Z' models. 
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where U = (u, c, t), D = (d, s, b), Vkm is the Kobayashi-Maskawa mixing matrix 
and Vr = (Vmj) = (V m * V^) with V R \j the unitary matrices which rotate the right handed 
quarks u Ri and dm from the weak eigenstate basis to the mass eigenstate basis. The model 
contains FCNC and their phenomenology has also been studied in detail, it is summarized 



in Ref. 17 



To have an acceptable neutrino mass spectrum and mixing, additional fields need to be 
introduced. For this purpose we introduce the Higgs triplet representations, A L : (1, 3, 1)(2) 
and Ar : (1,1,3)(2). For our model to be interesting for semileptonic B decays, it must 
have a light right-handed neutrino. This implies that the vevs and for the two Al and 
A R fields need to be small. In particular is not as large as in the usual Left-Right seesaw 
models. The neutrinos obtain Majorana masses from u A ' in addition to the Dirac masses 
induced by the non-zero vev of <fi. The resulting neutrino mass eigenstates (z/™, (z^ 3 ) c ) are 
related to the weak eigenstates (vl, u R3 ) by 



/ 




(9) 



Note that vl and have three components. 

Since only the third generation has a right-handed neutrino, Ul = (Uuj) is a 3 x 3 
matrix determined by the Majorana mass term for left-handed neutrinos. Url — {Urlis) 
and Ulr = (UiR3i) are 3x1 and 1x3 matrices, respectively, determined by the Dirac mass 
terms for neutrinos, and Ur = (Ur^) is a 1 x 1 matrix from the Majorana mass term for 
the right-handed neutrino. 

As already mentioned, we consider the scenario where the new right-handed neutrino is 
very light (even though it has a Majorana mass), as this is the only case relevant for B 
decays. If the third generation neutrino is mostly Dirac type, the elements in Ulr,rl can 
be large, of order one. In this case z/£ 3 is dominated by the weak eigenstate ur 3 . If the 
right-handed neutrino is Majorana, Ur can be large. If the right-handed neutrino has a 
mass negligible compared with the r mass, these two cases give the same results for our 
purposes. 

In general the charged leptons will also mix with each other, and their weak and mass 
eigenstates are related by: t^R = V[ R £™ R with £ = (e, /i,r) T . Here V[ R are the unitary 
matrices which diagonalize the charged lepton mass matrix. In the lepton mass eigenstate 
basis (dropping the superscript "m" on the fields), we have 

- ^(Ph-fUlajt* + PmT»Ui^jv)(mnfaW+ + cos^ w W'+) + h. c, (10) 

where 

U^ = UM, U& 4S = (Ui La Vl lj ), U e LRij = UumV^ , U e mj = U R33 V^. (11) 
U e is approximately the PMNS matrix. 
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In B meson decays the final neutrino flavor is not identified so it must be summed. For 
the processes involving left- and right-handed charged leptons, neglecting neutrino masses 
compared with the charged lepton masses, the final decay rates into a charged lepton £j are 
proportional to 

For t L3 : £ \U[ 3 \ 2 + \U RLj3 \ 2 = £ \W m Vt^ + \U RLl ,V[ l3 ? 
= (E Wlu\ 2 + lUnusfwL/ = E l^l 2 = 1 , 

i I 

For £ Rj : £ \Ui Rij \> + If/^f = E ItWl,/ + I^W^/ 



= (E l*W + |%3 3 | 2 )|^ 3 / = \V l mj ? ■ (12) 

i 

In obtaining these results we used the unitarity of U: Yj% Wlu\ 2 + Wrli^ 2 = 1> \Ulrzi\ 2 + 
1^33 1 2 = 1 and the unitarity of V[ R . 

All this implies that within this model, the leptonic and semileptonic decay modes of 
the B meson involving r-leptons will turn out to depend on the following combinations of 
parameters (two for each case q = u,c and for each lepton £): 

4 



\v £ 

\ v mi 


21 


V R qb 


2 






2 




^ = ^^ifel (l - (M) ( 1 + f^V IO' I ■ (13) 



The first term arises from direct W contributions and the second term appears through 
W — W mixing contributions. In the next section we summarize our results in terms of the 
combinations of Eq. [131 



III. LEPTONIC AND SEMILEPTONIC 5-MESON DECAY MODES INTO r- 
LEPTONS 

A. Inclusive semileptonic B decay 

The inclusive semileptonic decay rate is easily calculated to first order in the W — W 
mixing parameter using the free quark model. The differential distribution with respect 
to the dimensionless momentum transfer z = Mj y /m 2 , Mf v = (pb —p c ) 2 = (pe + Pu) 2 , can be 
written in terms of the combinations defined in Eq. [13] as 



dT{b^q£-u T ) Glm\ 



V qb \ 2 ^\{l,rlz) 



2\2 



W ' 2? 



dz 192vr 3 

• (z{r\ + r 2 q (z - 2) - 2z 2 + z + 1) - rj(z 2 + z(l + r 2 q ) - 2(1 - r 

_ 24 Fl- r - 

^± r Mix ' q 



(14) 
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in terms of the usual kinematic function 

X(x, y, z) = x 2 + y 2 + z 2 — 2{xy + xz + yz), 



(15) 



and of the ratios r q = m q /mb and 77 = m^jm^. For electrons and muons, where the lepton 
mass can be neglected, this can be easily integrated to obtain the total decay rate, 



r(6 -> ql'v) 



G\m\ 



V, 



qb\ 



m 



2\ 1 



Fw> F ( — § + - F^ ix G — § 
mi J 2 \mi 



1927T 3 

where F(r) and G(r) are the known kinematic functions given by 

F(r) = 1 - 8r + 8r 3 - r 4 - 12r 2 lnr 

G(r) = -8y/r (l + 9r- 9r 2 - r 3 + 6r(l + r) lnr) 



(16) 



(17) 



Notice that the term proportional to W — W mixing behaves as an anomalous cbW coupling. 
In the M W i » M w limit it corresponds to 



9nV c f 
9l Vcb 



(18) 



in the notation of Ref. [2l|. 2 

It follows from Eq. [TBI that, in the limit where W — W mixing is negligible, the CKM 
angles extracted from inclusive b decay are modified as follows 



V q f 



qb 



/ g R M w 



\9lM } 



\v e 


2 |^| 2 




v qb 


2 



(19) 



and would be different in decays into electrons or muons according to the mixing angle V^ 3e . 
This would need to be taken into account in a global fit using our model. Since there is 
no evidence of a difference in V q b as extracted from electron or muon modes we will take 
V^ 3e = Vr 3 „. Also, we are interested in new physics that affects mostly the r- lepton, so we 
will assume for now that V^ 3e = V^ 3 ~ 0, so that V^ 3t ~ 1. 

For decays into r-leptons, I = r, the kinematics is more complicated but Eq. [TH can be 
integrated to give 



T(b — > qr u T ) 



Gl ml 



1927T 3 

in terms of the kinematic form factors 



V r , 



qb\ 



Fl> F(r 2 q ,r 2 ) 



4fL G(r 2 ,r 2 T ) 



(20) 



F(x, y) = y/\(l- 7(x +y)- 7(x 2 + y 2 ) +x 3 + y 3 + xy(12 - 7(x + y)) 
+ 12 



x Log + y Log 



Ax 



Ay 



2 There is an error in Eq. 18 of Ref. 2l|. Instead of 2G it should say G/2 and this leads to —0.015 < 
Re(K«) < +0.C 
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— x 2 y 2 Log 



[1-x-y + y/Xf 
Axy 



G(x, y) = y/x [\/\ (l + x 2 + lOx - 5y{l + x) - 2y 2 ) 

{l-y)(l-y-VX)-x(l + y y 



+ 6x(x + (1 — y) 2 )Log 
— 6y 2 (l — x)Log 



x(l — x + y — VT) 

(l-x + y-y/X) 2 " 
4y 



(21) 



where we have used A as a shorthand notation for A(l, x, y) as defined in Eq. [15] The g B — > 
limit of Eq. [21] agrees with the known SM result [l8, 19]. Corrections to the free-quark decay 
model, of order 0(1/ m 2 ) have been calculated in the literature for the SM case [19| and can 
be added to our results for detailed comparisons with experiment. 

For our later estimates it will be useful to have the numerical results for the kinematic 
form factors. Taking m c = 1.25 and rrif, = 4.2 we obtain F = 0.069 and G = 0.014. 



B. Leptonic decay — > £ U£ 

The rate for the leptonic decay B^r — > (q = u, c) is given in terms of the 5 g -meson 

decay constant fg defined by 

< 0\qYl 5 b\B q (p) > = -if BqP » (22) 



as 



T{B± Put) = ^f 2 Bq m 2 m Bq (l-^f) \V qb \ 2 - 2 i&J (23) 

V m B q J 

For the leptonic decay only the axial- vector current contributes. In the limit of no W — W 
mixing, Eq. [23] reflects that there is no interference between the W and W contributions 
due to the different Dirac structure in the lepton leg. 



C. Semileptonic decay B — > Dtv 



To calculate the exclusive semileptonic decay rate for B — > Dtv we use the standard 
parameterization for the form factors with q = p — p' 



< D( P ')\cY b\B( P ) > = U(q 2 )(p + p'r + f-(q 2 )(p-p'Y- 



(24) 



In terms of these form factors and the kinematic function of Eq. [15] the differential decay 
rate can be written as 



dT{B Dtv) G 2 m\ yJX(l, r%, z){z - r 2 ) 2 



dz 



192vr 3 

f2 fn„.2 2\ 



2z 3 



fi(2X(l } r 2 D} z)z + 



l,r 2 z) + 3(l-r 2 



D, 



/+/_(6(r* - \)r\z 



I Kb 1 2 ( F W> + 2 ^Mix) 



(25) 
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where z = Mj v /M B = q 2 /M B) r D = M D /M B and r e was defined above. 

It is interesting that this expression, Eq. |25l is proportional to the SM result. This occurs 
because the decay is only sensitive to the vector form factor which has the same sign for 
both W and W contributions. This result allows us to write the decay rate in this case as 

Y{B Dtp) = T(B -> Drv) SM + 2 F Mix ) (26) 

and no further knowledge of the form factors Eq. |2U is needed in order to compare with the 
experimental results. 



D. Semileptonic decay B — > D*tv 

To calculate the exclusive semileptonic decay rate for B — > D*tu we use one of the 
standard form factor parameterizations 20|, modified so that all the form factors are di- 
mensionless: 

<D*(p')\c^ l5 b\B(p)> = -iM B f(z)e^-i^-(a + (z) (p + p>y + a_(z) (p - p'Y) 

Mb 

< D*(p>)\cY b\B(p) > = 9 M e ^ el (p + pXip-p')?. (27) 

Mb 

In terms of the dimensionless variables z = Mf u /M B = q 2 /M B , tb>* = Mb,*/ Mb and = 
Mij Mb we obtain 



dT(B — >■ D*tu) G\m\ 

cTz ~ 1927T 3 



B 



\V qb \ 2 [Fyy, (F v + F A ) + 2 F« ix (F v - F A )} (28) 



where the two kinematic form factors are now given by 

A 3 / 2 

Fv = 9{zf — (2z + r 2 ) (r 2 - zf 



r^z- 



2\z[X+ l2r 2 D *z) + 2r 2 (\ + 3r 2 D *z)) f(z) 2 

+ 3Xr 2 z 2 a„(z) 2 + A (2\z + r] (A + 3(r£* - l) 2 )) a + (z) 2 
+ 6\r 2 z((r 2 D * - 1) a+(z) a-(z) - f(z) a_(z)J 

- 2A (2z (r 2 D , +z -l)+r 2 (ir 2 D , + z - 4)) f(z) a + {z)\ (29) 

where A in this expression is short for A(l, rf^*, z). This expression indicates that for this 
mode the new physics is not simply proportional to the SM because the contributions from 
the vector and axial-vector form factors enter with different signs in the term arising from 
W-W mixing. 

To judge the relative importance of the two terms: W contribution vs W — W mixing 
contribution, we provide a rough numerical estimate for their respective contributions to the 
decay rate: 

d 2 m 5 

T(B -> DVi/) = \V qb \ 2 (0.062 F*,, - 0.11 F« ix ) (30) 

To obtain this estimate we converted the form factors to the usual HQET parameterization 



201 ] and used = 1 for the Isgur-Wise function over the full w = vb • Vr>* range. 
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IV. NUMERICAL RESULTS AND CONCLUSIONS 



The data on semileptonic B decay into r leptons Eqs. [U [2] can be fit to our model in 
terms of the two combinations of constants defined in Eq. [T3J For example, if we use the 
central value of each experiment we obtain 

BaBar => Ffc, = 1.40, i^ ix = 0.043 
Belle Fyy> = 1-46, i^ ix = -0.14 (31) 

If instead we use the more precise BaBar result at the one sigma level we find 

1.2 < Ffr, <1.6 
0.009 < F^ ix < 0.08 (32) 

To see whether these numbers can be produced in our model let us recall the existing 
constraints. The model described earlier is geared for an enhanced ratio of right-handed to 
left-handed couplings gn/gi- The constraints that LEP2 data impose on this ratio can be 
roughly summarized by the relation [l5| 

* £ ^, (33) 
9l M w 

whereas the requirement that the new gauge interaction (SU(2) R ) remains perturbative 
implies that g R ^ 10 gi- 

The FCNC generated in this model have also been studied in detail and summarized 



in Ref. [171] . There, an ansatz was produced to understand the smallness of the observed 
FCNC, and this ansatz gives us specific predictions for the quark-mixing angles that appear 
in Eq. [131 Specifically, we found that 



which then lead to 



Vhk ~ fa, Vmc ~ V cb , V£ tu ~ V ub (34) 



V Rub ~ V u b V Rcb ~ V cb . (35) 



Finally, as mentioned earlier, we will assume that the charged lepton mixing matrix is nearly 
diagonal so that V^ 3r ~ 1. Under these conditions the fit of Eq. [321 requires 



n 7 < 9rM w < n n 
U "' ~ 9lM w , ~ U ' J 

0.005 < ^ja <o.06 (36) 



The first of these conditions can be satisfied within the constraint Eq. [33]and the requirement 
that the new interaction remains perturbative provided My/ 1 ~ 1 TeV. This is a range that 
should be accessible to LHC, providing new motivation to search for a non-universal W . 
The second condition requires that we re-examine bounds on W — W mixing, in particular 



the combination £ e // = £w9r/9l- This combination is constrained by b — > s'j [15|, |22) and 



we applied this constraint to our model in ref. [15(. We can update that result using a recent 
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HFAG average B(b — > X s '-f) = (3.55 ± 0.25) x 10 4 23[ for the experimental number and a 
recent NNLL SM calculation that gives B(b -> X s ^) = (3.15 ± 0.23) x 10" 4 [24J. Assuming 
that the new physics interferes constructively with the SM, at the 3a level the allowed range 
becomes 

- 1.3 x 10~ 3 < £ e // £ 2.7 x 10~ 3 (37) 

For destructive interference between the new physics and the SM there is a second, narrow, 
allowed range near £ e ff = —0.028. We see that the constraint from b — > 37 limits the effect 
of mixing to be somewhat lower than what is required for the fit Eq. [321 at la. 

Our conclusion from this study is that our model can accommodate an enhancement in 
the ratios R(D) and R(D*) of the order reported by BaBar and Belle (through Fyy>), but 
that it predicts approximately the same enhancement for both ratios once we incorporate 
the constraint from b — > 57 on F^. 

Our model also leads to predictions for related modes. For the leptonic B — >■ tu decay 
we find 

T(B~ -> r~u) 



T(B ->■ r u T ) S M 



F W' - 2 ^Mix- ( 38 ) 



With the constraints discussed above, this implies that a number as large as implied by the 
BaBar data in Eq. [4] is possible. However, the mixing angles involved are different, in this 
case Vjiub instead of Vr C 6, so this enhancement is not necessarily equal to the one occurring 
for R(D) and R(D*) although it is expected to be of the same order. 
The inclusive semileptonic decay b — > cr~v can also be enhanced 

r(6m) 0.8*5- (39) 



r(6 -> CT U) SM 

and the constraints described above from 6—7-57 restrict this factor to be approximately 
the same as the one occurring in the exclusive modes. 

Finally, the model also leads to a correlation between the enhancement observed in 
R(D, D*) and the leptonic decay of B~ 

T(B- -> t-v) 



T(B- T -v T ) 



SM 



Fw> ~ 2 F^. (40) 
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